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ABSTRACT 

Recent observations of coronal- loop waves by TRACE and within the corona as a whole by 
CoMP clearly indicate that the dominant oscillation period is 5 minutes, thus implicating the 
solar p modes as a possible source. We investigate the generation of tube waves within the 
solar convection zone by the buffeting of p modes. The tube waves — in the form of longitudinal 
sausage waves and transverse kink waves — are generated on the many magnetic fibrils that lace 
the convection zone and pierce the solar photosphere. Once generated by p-mode forcing, the 
tube waves freely propagate up and down the tubes, since the tubes act like light fibers and form 
a waveguide for these magnetosonic waves. Those waves that propagate upward pass through 
the photosphere and enter the upper atmosphere where they can be measured as loop oscillations 
and other forms of propagating coronal waves. We treat the magnetic fibrils as vertically aligned, 
thin flux tubes and compute the energy flux of tube waves that can generated and driven into 
the upper atmosphere. We find that a flux in excess of 10 5 ergs cm" 2 s" 1 can be produced, 
easily supplying enough wave energy to explain the observations. Furthermore, we compute the 
associated damping rate of the driving p modes and find that the damping is significant compared 
to observed line widths only for the lowest order p modes. 

Subject headings: MHD — Sun: corona — Sun: helioseismology — Sun: magnetic fields — Sun: 
oscillations 



Introduction 



It is now well established through observations by the Solar and Heliospheric Observatory (SOHO), 
the Transition Region and Coronal Explorer ( TRA CE) , and Hinode that the solar atmosphere supports 
the propagation of a variety of magnetohydrodyn amic (MHD) waves. Slow magnetosonic waves have been 



implicated in the observat ions of coronal plumes ( Ofman et al. 1997 : Ofman. Nakariakov fc DeForest 19991: 



DeForest fc Gurmanlll99ct) and in corona l-l oop oscillations (e.g ., Berghmans fc Clette 19991 : be Moortel. Ireland fc Walshl 
2000': be Moortel. Hood fc Ireland! Eooj ). IWang et all \20oA) interpreted oscillatory behavior seen in very 



1 This paper was originally published in the Astrophysical Journal on 10 April 2008, vol. 677, pp. 769-780 
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hot (> 5 MK) coronal loops by SO-ffO/SUMER, as standing slow magnetosonic waves of long wavelength. 
A propagatin g intensity oscillat ion with 5 minute periodicity in a "nonsunspot" coronal loop has also been 
reported by Marsh et al. ( 20031 ). They estimate a speed of 50-195 km s _1 for this wave and suggest, on 
the basis of Coronal Diagnostic Spectrometer (SOffO/CDS) data, that it is also present at chromospheric, 
transition region, and coronal temperatures. 

Fast magnetosonic waves ha ve also been invoked to explain coronal and chromospheric oscillation ob- 
servations. Williams et al. ( 20021 ) report intensity fluctuations with a period of 6 s along and within coronal 



loops. The pr opagation speed is estimated to be 2 x 1 3 km s x , suggesting the possibility of fast magne- 



tosonic waves. 



Verwichte. Nakariakov fe Cooper ( 20051 ) have recently claimed to have se en propagating fast 



kink waves in a postflare supr a-arcade. Transverse loop oscillations observed by TRACE (jAschwanden et al 



1999; 



Nakariakov et al. 19991 ) have also been lin ked to fast kink waves. Finally, there have been sugges- 



tions that the large-s cale EIT waves discovered by Thompson et al. ( 19981) a r e the result of propagating fast 
magnetosonic waves (|Wang||2000l ; iMurawski. Nakariakov fc Pelinovsky||200ll ; lOfman fc Thompsonll2002l ). 

Recent observations by the Co ronal Multichannel Polarime ter (CoMP) instrument (jTomczyk et al. Il2007h 
and TRACE (see the overview bv lDe Moortel fc Rosnerll2007l) clearly indicate that coronal oscillations are 
ubiquitous and have amplitudes that peak at the p-mode frequency band. In the high corona, the CoMP 
instrument detects strong fluctuations of the line-of-sight Doppler velocity with a broad peak at 3.5 mHz. 
Correspondingly weak intensity oscillations are seen at the s ame f requency. This indicates that the observed 
waves are nearly incompressive, and as such lTomczvk et al.l (|20071 ) identified the waves as Alfvenic. However, 
in a low-/3 plasma, fast magnetosonic tube waves or kink waves are also nearly incompressive. Therefore, we 
posit that these observations in the high corona may also be fast kink waves. 

On the other hand, the loop oscillations observed by TRACE have been identified as slow w aves (e.g., 
Berghmans fc Clettel[l999[ be Moortel. Ireland fc Walshlboool : be Moortel. Hood fc Irelanj2002h . The os- 
cillations associated with sunspot loops have peak power at a period of 3 minutes (or a frequency of 5.5 
mHz), while loops not associated with sunspots have peak power at a period of 5 minutes (or 3.3 mHz; 
De Moortel et al.l 120021 ; ICenteno. Collados fc Truiillo Buenol 12006). Therefore, both slow and fast waves 
with frequencies associated with p-mode oscillations have been observed. The p modes, and possibly solar 
granulation, are clearly implicated as the ultimate source of many coronal fluctuations. In this paper, we 
concentrate on p modes as the driving mechanism. 

If p modes are in fact responsible, the p modes must somehow manage to leak into the upper atmo- 
sphere, despite the fact that they become evanescent above the photosphere. Some have suggested that 
the leakage is accomplished by a "ramp effe ct," by which an inclined magnetic field effectively reduces the 



acoustic cutoff frequency of the atmosphe re ( De Pontieu. Erdelyi fc De Moortell 2005 : Hansteen et al. 20061 : 



Mcintosh fc Jefferies 2006 : Jefferies et al. 20061 ). However, we point out that the cutoff frequency for both 
the slow sausage wave and the fast kink wave can fall be l ow th e nonmagnetic acoustic cutoff frequency 
(see Roberts fc Webb 1978| and Musielak fc Ulmschneider 2001 1 for a derivati on of the cutoff for sausage 
and kink waves, respectively). Therefore, a flux tube can act as a waveguide ( Robert J 1981 ). permitting 
efficient penetration of the acoustic barrier presented by the photosphere. It may be possible that vertically 
aligned flux tubes are capable of carrying a sufficient flux of tube waves into the atmosphere to explain the 
observations. 

We suggest the following mechanism. Within the solar convection zone, p-mode oscillations buffet the 
many magnetic fibrils that thread through the photosphere into the chromosphere and corona. This buffeting 
excites both sausage and kink waves on the fibrils, which then propagate both up and down the vertical field 
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lines. Those waves that travel downward are lost in the convection zone, whereas those that propagate upward 
pass through the photosphere into the upper atmosphere, where they manifest as coronal-loop oscillations 
and upward-propagating waves. 

In this paper we estimate, by a semianalytic method, the atmospheric energy flux of sausage and kink 
waves that can be generated by this mechanism. Our calculation is a direct extension of the work by 
Bogdan et al. ( 19961) . hereafter referred to as BHCC96. We demonstrate that / and p modes efficiently 



generate tube waves and that an energy flux in excess of 10 5 ergs cm~ 2 s _1 can be driven upward through 
photospheric levels. Whether, such a flux can survive passage through the chromosphere and transition 
region into the corona is left to subsequent work. 

The paper is organized as follows. In §2 we describe the stationary configuration of the flux tubes and 
the equilibrium atmosphere in which they are embedded. §3 details the p-mode oscillations, and §4 derives 
the driven oscillations of the thin flux tubes. In §5 we present our primary results, the energy flux of tube 
waves and the damping rates of the p modes that those fluxes engender. §6 provides a discussion of our 
findings, where we compare our derived atmospheric energy fluxes with those obtained by observations. 



2. The Equilibrium 

We examine the excitation of waves on a thin magnetic fibril through the buffeting of the tube by 
acoustic oscillations within the solar convection zone. The fibril threads through a field-free atmosphere, 
and its axis is aligned with the force of gravity. We assume that in the absence of the acoustic wave field, the 
fibril and surrounding atmosphere are static and in equilibrium. In the following two subsections we describe 
the field-free atmosphere in which the fibril is embedded and the equilibrium configuration and structure of 
the fibril itself. 



2.1. The Atmosphere 

We restrict our attention to acoustic waves with wavelengths much shorter than the radius of the 
Sun; therefore, we may ignore the curvature of the solar surface. This allows us to make the simplifying 
assumption that the atmosphere is plane-parallel with constant gravity pointing downward, g = —gz. The 
height coordinate z is defined to increase upward. We model the solar convection zone with an isentropically 
stratified atmosphere that is excised above a height z = —zq, known as the truncation height. Below this 
height, z < —zo, the atmosphere is polytropic and the gas pressure, mass density, and sound speed are 
power-law functions of height, 



/ \ m - 1 / \ m+1 

m+1 \ z J \ z J 

d{z) = - 9 ± . (2.3) 
m 

In the preceding equations, the subscript "e" is used to denote the nonmagnetic atmosphere that is "external" 
to the magnetic fibril. The constants po and P Q are the surface values of the mass density and gas pressure. 
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The constant m is the polytropic index, which has a special value that depends on the ratio of specific heats 
7, because the atmosphere is isentropically stratified, m = 1/(7— 1). 

An extremely diffuse and hot gas exists above the truncation height. For this work, we adopt the limiting 
case of a vacuum (p — > 0) with infinite temperature (T — * 00), such that the pressure remains finite. At the 
surface, where z — — z , the pressure of this "hot corona" matches the pressure of the lower atmosphere (the 
model convection zone). 

Following BHCC96, we specify the truncation depth zq and the surface mass density po (and therefore 
the surface pr essure Pp) by matc hing the surface layer of our model, z = —zq, with the T5000 = f level of a 
solar model bv lMaltbv et al.l (|l986h . At this layer in the solar model g = 2.775 x f 4 cm s -1 , p = 2.78 x 1 0~ 7 



g cm~ 3 , and Po = 1-21 x 10 5 g cm -1 s _1 . We adopt m = 1 .5 (appropriate for 7 = 5/3), which is a reasonable 
approximation to the bulk of the convection zone. With this choice of polytropic index, the truncation depth 
and surface sound speed become 392 km and 8.52 km s _1 , respectively. 



2.2. The Magnetic Fibril 

For simplicity, we assume that the magnetic fibril threading the truncated polytrope obeys the following 
properties: 

f . The flux tube is straight, untwisted and axisymmetric. 

2. The tube has a circular cross section. 

3. The tube is vertically aligned and coincident with the z-axis. 

4. The tube has a potential magnetic field lacking internal currents. 

5. The boundary of the flux tube is marked by a current sheet. 

6. The tube is thin, meaning that its radius is significantly smaller than any scale length in the atmosphere, 
including the density scale height and the wavelength of incident acoustic waves. 



The last of these assumptions bears some further discussion. A sufficiently thin flux tube is unable to 
support internal forces and structures beyond hydrostatic balance. Thus, the gas pressure imposed by the 
external medium at the flux tube boundary must be matched by a total pressure that is constant with radius 
across the tube. This same argument holds for temperature as well. The temperature inside and outside 
must be identical (as are the pressure scale heights), because any initial difference would be destroyed by 
thermal and radiative diffusion. Continuity of total pressure across the flux tube interface requires that the 
magnetic pressure also has the same scale height. Hence, the value of the plasma parameter j3 — 8irP/B 2 
must be constant with height within the tube. 

This simple statement of constant (3 within the tube is a strong constraint on the geometry of the 
flux tube. Since the gas pressure in the solar atmosphere decreases rapidly with height, the magnetic field 
strength within the tube also decreases rapidly, and the flux tube must flare dramatically near the surface. 
Thus, at some height in the atmosphere, the approximation that the tube is thin must break down. The 
magnetic field outside of sunspots is generally found in the form of small magnetic elements. For such 
elements, the flaring that occurs with height becomes important within the chromosphere where neighboring 
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flux elements begin to collide, thereby forming a magnetic canopy. The thin-flux-tube approximation cannot 
be applied consistently within the upper atmosphere. In this work, we circumvent this problem by truncating 
the atmosphere and only studying wave propagation in the region where the flux tube remains thin. As will 
be discussed more fully in §4.3, the details of the tube geometry and thermodynamic structure within the 
upper atmosphere can be mimicked by an appropriate choice of boundary condition at the surface, z = — zq. 

For thin flux tubes of the sort described above, radial variation of the field strength can be ignored and 
the flux tube is completely defined by its total magnetic flux O and the plasma (3. The tube's internal gas 
pressure P(z), mass density p(z), axial field strength B(z), and cross-sectional area A(z) are given by the 
equations 



P(z) = ^fi^eW, (2-4) 
P{z) = ^-yPc(z), (2.5) 

For a flux tube with (3 = 1.0 embedded in a m = 1.5 polytrope with a surface mass density of po = 2.78 x 10~ 7 
g cm~ 3 , the magnetic field strength at the surface, z = —zq, is 1.2 kG. 



3. The p-Mode Oscillations 

The nonmagnetic atmosphere surrounding the magnetic fibril supports acoustic oscillations that arc 
trapped in a wave guide just below the surface. The upper boundary reflects upward-propagating waves and 
the increasing sound speed with depth refracts downward-propagating waves back towards the surface. Since 
the atmosphere is isentropic, the acoustic modes of this waveguide — the / and p modes — have a particularly 
simple form and can be expressed using a displacement potential $, 



5P C 

8p e 



<9 2 $ 
Pc d 2 <$> 

a dt 2 ' 



(3.1) 
(3-2) 

(3.3) 



The vector £ e is the fluid displacement, 8P C is the perturbed gas pressure, 8p c is the perturbed mass 
density and c e is the sound speed. Since the tube is thin, it is incapable of supporting a temperature 
differential with its surrounding; thus, the sound speed inside and outside the tube are the same function of 
height c(z) = c c {z). We drop the subscript on c c from here on. 
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3.1. Wave Equation 

Equations (|3.1[) - (|3.3[) can be combined into a single partial differential equation for the displacement 
potential, 



Equation (|3.4[) supports plane- wave solutions of the form 



(3.4) 



$(x,t) = Ae~ lu ' t e lkx Q{z) 1 



(3.5) 



where A is the complex wave amplitude, u is the temporal frequency, k is the wavenumber, and Q(z) is the 
vertical eigenfunction; Q(z) is dimensionless and A has units of length squared. We have assumed without 
loss of generality that the acoustic wave propagates in only one horizontal direction, the x-direction. Direct 
substitution into equation (|3.4j) produces an ODE for the eigenfunction Q[z), 



^-»i^-*v»e (2 ) = o. 



(3.6) 



This equation can be cast in a useful dimensionless form by the substitutions 



2 _ mu) Zq 
9 



w — —2kz, X = 2kzQ 1 



X ' 



resulting in the equation 



dw 2 



m d 
w dw 



Q(w) = 0. 



(3.7) 



The two solutions to equation (|3.7p involve Whittaker's functions W K ^(w) and M K ^(w) fsee lAbramowitz fc Stegun 



19641) 



lM K(l (w) 



(3.8) 



where fi = (m + l)/2. Whittaker's W function vanishes as w — > oo, but is poorly behaved at the origin 
w — 0, whereas Whittaker's M function is well behaved at the origin and diverges exponentially asm-* oo. 



3.2. p-Mode Boundary Conditions 

As boundary conditions we require that the solution vanishes deep in the atmosphere as w — > oo, and we 
require that the Lagrangian pressure perturbation vanishes at the free upper surface z = — z (or equivalently, 
w = X = 2kz ). The first of these requirements makes Whittaker's M function unsuitable. The second of 
these boundary conditions places a restriction on the allowed values of k (or equivalently A). Therefore, for a 
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fixed value of the frequency, the eigenvalue quantitized and satisfies the following set of equations, 

which simply restate that the Lagrangian pressure perturbation vanishes at the surface, 



d_ 

dw 



Qn(Xn) = 0, 



(3.9) 



where, 



(3.10) 



After using recursion relations for Whittaker's W function fsee lAbramowitz fe Stegunlll964l ). the last two 
equations can be reduced to 



W Kn ^ +1 (X n ) =0. 



(3.H) 



The lowest order solution with n = is the fundamental mode or / mode, whereas the higher overtones 
n > correspond to the p modes. 



3.3. Acoustic Mode Energy 



The energy density of an acoustic wave propagating through an isentropic media is given bv lBrav fe Loughhead 



(119741) 



15P? 



2 Pe C 2 



(3.12) 



where the vector £ e is the partial derivative of £ with respect to time. 



After taking temporal averages of equation (|3.12p . integrating over height, and multiplying by the surface 
area of the Sun, one finds the following expression for the energy contained in /- and p-mode oscillations, 



E„ 



Mi 



4^ 2 ^^-^Mi 
4m A™ zl 



dw 



( dQ n (w) 
\ dw 



W 



(3.13) 
(3.14) 



Excitation of Tube Waves 



The magnetic fibril is buffeted and driven by /-mode and p-mode waves within th e convection zone . 
Magnetosonic oscillations of thin flux tubes have been ex amined in detail by many authors ( Spruit][l981 . 1984 : 
Stix 1991 : Rvutova fc Priest 1993t Bogdan et al. 1996 ). and can be fully described by two types of waves: 
kink and sausage. Sausage waves are axisymmetric pressure pulses that produce displacements that are 
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primarily parallel to the magnetic field, £u(z,t). Kink waves produce perpendicular displacements, £j_(z,t), 
and magnetic tension and buoyancy are the primary restoring forces. 



4.1. Thin-Flux- Tube Equations 

The sausage waves and the kink waves are driven respectively by the overpressure and the transverse 
velocity imposed on the outer surface of the flux tube by incident acoustic waves. Using the formulation of 
BHCC96, the sausage and kink waves can be described by the equations 



dt 2 T dz 2 2c 2 9zJ qi pVldzdt 2 ' K ' } 

2 2 d 2 , jgcid 1 _ o c c 2 K 



dt 2 K dz 2 2c 2 dzj^ pV 2 dxdt 2 ' [ ' 

where Va is the Alfven speed, ct is the cusp or tube speed, and ck is the kink speed. This later speed is also 
referred to as the "mean" Alfven speed, since it is the density-weighted mean of the Alfven speeds inside 
and outside the tube: 

2 B 2 2 c 2 V^ 2 B 2 



Vt — -— r z — A r z 



4^p' ^ c 2 + ^' 47r(p + p e )' 

Equations (|4.ip and (14. 2p describe the forced oscillations of sausage and kink waves on a slender tube. 
Clearly, in the limit of low 0, the saus age waves are simply slow magnetosonic waves propagating along the 



thin tube (see lRoberts fc Webbl [1978j for a derivation in an atmosphere with generic stratification). In this 
same limit, the kink waves are fast magnetosonic waves, with an enhanced effective density due to the fact 
that transverse motions of the tube must push and pull external fluid. However, in the parameter regime 
that corresponds to photospheric magnetic elements, the plasma is roughly unity and the distinction of 
fast versus slow magnetosonic waves proves less useful. In any event, the sausage wave is primarily a pressure 
wave, whereas the kink wave is principally a tension wave. 

The shaking of the fibril by p-mode oscillations appears as the forcing on the right-hand sides of equa- 
tions (|4.ip and (|4.2[) where the pressure perturbation and horizontal displacement have been replaced by 
appropriate derivatives of the displacement potential. The forcing for the sausage waves arises from the 
requirement that the total pressure is continuous across the tube boundary. The kink waves are driven by 
horizontal motions of the external media. 

These equations have similar form and can be written compactly using the notation, 

-{icr = U{S), (4.3) 



ds 2 s ds s J zq 

where s — — z/zq is a dimensionless depth and a represents either || or _L, corresponding to the sausage and 
kink waves, respectively. The quantity e a takes on the two values, 

_ 2m + f3(m + 1) _ (m+ l)(2/3+ 1) 

2m 2m 
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The forcing function / CT (s) depends on whether we are driving kink or sausage waves, 



_ (m+l)(/3 + l) ^ dQ n (s) 

; " (s) = 2^ T - ' (44) 

.(to+IX/3 + 1) A„z/ 2 

/x(s) = i Q n (s). (4.5) 

2m s 

The homogeneous solutions to equation (|4.3[) can be expressed in terms of Hankel functions, 



Vv(s) 



(4.6) 



Us) 



(4.7) 



For the chosen time dependence (see eq. |3.5j ). the tp a solution represents a downward-propagating wave, 
while the a solution is an upward-propagating wave. Note that the two solutions are simply complex 
conjugates of each other, thus allowing all further equations to be expressed using solely i/j a . 

The solution to the driven problem is constructed from the homogeneous solutions using a Green's 
function formulation, 



Z Zq 



fin 



dr r" +1 ^*{r)U{r) 



+ 6 a (s) drr^ +1 ip a {r)f a (r) 



(4.8) 



This solution was constructed to satisfy specific boundary conditions. Deep in the atmosphere, as 
s — ► oo (or z — > — oo), £ CT becomes proportional to i/v, thereby satisfying a causal radiation condition. At 
the surface s = 1 (or z = ~zq), the solution is a mixture of upward- and downward-propagating waves. The 
exact proportion of ij} a to 9 a in this mixture is completely specified by the boundary-condition parameter 
fio- . The boundary condition represents both the reflection that occurs directly at the upper surface and the 
combined effect of the entire upper atmosphere z > — zq, which is not included explicitly in the calculation. 
Any physical choice of boundary condition can be satisfied by an appropriate value of fio-. For example, 
fi ff = enforces a radiation condition, where only the upward-propagating wave is present at the surface. 



The integrals that appear in equation (|4.8p form the basis for all of the derived quantities that appear 
later. The p-mode eigenfunctions Q n are purely real functions, and therefore the driving functions f a are 
either purely real (for the sausage) or purely imaginary (for the kink) . This allows us to express both integrals 
compactly for both the sausage mode and the kink mode, 



Ms) = J drr^ +1 Mr)U(r), 
T a = lim J a (s) = / drr"+V CT (r)/ CT (r), 

dr r» +1 i> a (r)fa(r). 



(4.9) 
(4.10) 
(4.11) 
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We call X a the interaction integral between the p mode and the respective tube wave. Using these 
definitions explicitly, the displacements for the two types of waves take on the forms 



£y(«) 

Note the difference in sign in front of J* in the first term in the curly braces. This difference arises because 
/|| is purely real and f± is purely imaginary. 

In equations (|4.12p and (|4.13p the term ipa-^T* represents the upward-propagating wave generated by 
the driver, while the term ip* [2 a — J a \ represents the downward wave directly generated by the driver. The 
term t/> CT Q CT is a downward-propagating wave arising from reflection off the upper surface, plus any waves 
propagating downward from the upper atmosphere through the upper surface (perhaps caused by a reflection 
in the chromosphere or corona). 

Note that the energy flux of waves generated by the driver and propagating away from the driving region 
is the same both upward and downward. This can be seen by examining only those terms in equations (|4.12|) 
and (|4.13[) that represent wave components directly generated by the driver. The downward component is 
evaluated as s — > oo, and the upward component at s = 1, 



iir A n 
'~2~Zo 



(a) 0|| + ^( s ) +^( S ) in-Jj|(s) 



l ^U x { S ) <> L -J ,.s) r (.s) I ■■:/(,) 

l Zq 



(4.12) 
(4.13) 



^'^^W, (4-14) 

Z Zq 

e p) = ~—r^. (4.i5) 

2 z 



This is a useful fact that will be exploited later. 



4.2. Energy Flux of Tube Waves 



At any point along the tube, the energy flux is given by the following expression (jBrav k, Loughhead 



1974), 



'-(" + ^)*-^ fl - 

The perturbed magnetic field for sausage waves is given by equation (A8) in Appendix A of BHCC96. The 
perturbed field for the kink waves can be derived directly from the induction equation of MHD, 

— = 2TTP SPc + ^TT^^" 9 { p ~^TTP Pc ) Ch (4 - 17) 

SB X = B^. (4.18) 

az 
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Direct substitution of equations (|4. 1 2[) . (|4.13[) . (14. 1 T[) and (|4. 18|) into equation (|4. 16|) produces the 
following vertical energy fluxes at the two boundaries. As s — > oo, 



and at s = 1, 



where, 



Fn 



F± 



7/3 ngpo^ 



2 + 7 /3 4(m+l)(/3 + l) 



4(m + 1)09 + 1) zl 



Qii +Xn 



F = 7/3 Trgppcj \A n \ /, ,2 

11 2 + 7/3 4(m+l)(/3 + l) IJ III 



4(m +!)(/?+!) z£ 



fill +<S , 



(4.19) 
(4.20) 



(4.21) 
(4.22) 



Oil 



<>;, 



(4.23) 



Qn = Qn(s = 1) = Q„(w = An). 



(4.24) 



Note that the flux passing both down and up the tube are functions of the boundary condition that is 
applied at the surface s = 1 through the fl a parameter. Furthermore, the sausage mode also possesses a 
term S in the energy flux at the upper surface s = 1 which arises because the driver, equation (|4.4[) . does not 
vanish at the surface. This term is essentially the p-mode eigenfunction multiplied by the real part of the 
vertical displacement within the tube. Depending on the phasing between these two quantities, this surface 
driving term can take on positive or negative values. 

Furthermore, the flux deep in the atmosphere s — » oo increases with depth F ~ s M+1 , but it does so at 
the same rate that the cross-sectional area of the tube decreases, 



Ms) = A™+D(/?+i)y /2 es -(, +1 ) 

V &ngz po J 



(4.25) 



Therefore, below those layers in which p modes drive waves, the rate at which energy passes down the tube 
is constant with depth, 



Eu 



7/3 TrgpoioAs \A n \ 



2 + 7/ 3 4(m+l)(/3 + l) zl 



Sill + In 



(4.26) 
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E, = 



irgp ujA s \An\' 



4(m+l)(/3 + l) 



± 



(4.27) 



In these two equations, E is negative for energy escaping the p-mode cavity, or in other words, if the energy 
flux is downward. 

At the upper surface s = 1 a similar result is obtained, 



E\\ — — 



E± 



id 



irgp ujA s \A n Y 



2 + 7/?4(m+l)(/3 + l) zl 

Trgp uiA s \A n \' 
4(m +!)(/?+!) zl 







■5 , 



(4.28) 
(4.29) 



In these expressions, E is once again negative for energy escaping the p-mode cavity. However, in this case 
this requires that the flux is upward into the upper atmosphere above our model convection zone. 



4.3. Boundary Conditions 

We apply two different boundary conditions at the surface and track the results of each. The first 
boundary condition is the requirement that the stress vanishes at the upper surface s = 1. Of course, for 
this boundary condition, the energy flux through the upper surface must be identically zero. We include this 
boundary condition because it is the same boundary condition adopted by BHCC96. The second boundary 
condition that we apply maximizes the energy flux that passes up through the upper surface. This boundary 
condition allows us to compute an upper limit on the amount of tube-wave energy that is driven into the 
upper atmosphere. 



4-3.1. Stress-Free Condition 

For the sausage mode, the stress-free condition is enforced by requiring that the divergence of the 
displacement vector vanishes at the upper surface. This is equivalent to setting the Lagrangian pressure 
perturbation to zero. BHCC96 have shown that for the sausage waves, 

After inserting 5P e and and after substantial manipulation, one finds that V • £ = requires, 

n„ = j^±m±A^.3. Xh ( 4.3i) 

«ii = ^Kii(. 2v VH) + ^+ 1 )^+ 1 ) H l 1 H 2 ^)- ( 4 - 32 ) 
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Fig. 1. — Energy contained in solar p modes with harmonic degr ees in the range 3 < £ < 150 and orders n > 



Diamonds indicate the measurements of lKomm. Howe fc H ill (2003). As can be clearly seen, the energy is primarily a 



function of frequency. The solid red curve is a maximum-likelihood fit to the data using the functional form appearing 
in eq. (|5.3[1 . 



For the kink waves, BHCC96 imposed the boundary condition at s = 1 that the magnetic tension force 
vanishes or equivalently, that the second derivative of the displacement with respect to height vanishes. We 
do not follow their example, for the simple reason that the horizontal force equation for the kink mode 
possesses two other terms, the buoyancy force and the external p-mode forcing. In equation (|4.2p the term 
with the second derivative with respect to height arises from the magnetic tension, while the term on the 
left-hand side with the single derivative with respect to height is the buoyancy force. The external forcing is 
the right-hand side of the equation. Neither the buoyancy nor external forcing vanish when the tension force 
is zero. However, setting the displacement itself to zero at the surface is sufficient to specify a net horizontal 
force of zero. Therefore, this is the proper stress-free boundary condition that we apply. Evaluation of 
equation (|4. 13|) at s = 1 reveals that this stress-free boundary condition requires 

= - m V ' Ji. 4.33 

Direct substitution of the values of fl a given by equations (|4.31[) and (|4.33[) into the relevant energy- 
transfer rates, equations (|4.28|) and (|4. 29[> . reveal that the energy flux through the upper surface is zero, as 
expected. 



4-3.2. Maximal-Flux Condition 

The maximal-flux boundary condition allows us to place an upper limit on the energy contained within 
coronal-loop waves that are generated by p-mode forcing in the solar convection zone. We specify the value 
of £1 CT by maximizing the energy flux shown in equations (|4.2ip and (|4.22p . 
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Fig. 2. — Energy flux of tube waves driven up through the model photosphere. The solar surface is peppered with 
identical thin flux tubes with j3 = 1.0 and a combined filling factor /. The fluxes shown are divided by this filling 
factor. The different curves correspond to the flux driven by different order p modes: / (black), p\ (red), P2 (green), 
P3 (dark blue), pi (yellow), ps (aqua), P6 (violet), and P7 (ora n ge). T he range of frequencies for which the curves are 
solid corresponds to the window where iKomm. Howe fc Hill measured p-mode energies. Frequencies where 

the curves are dashed indicate extrapolations of the measured p-mode energy. These curves are for the case where a 
maximal-flux boundary condition has been applied at the photosphere. 



"« " J "' + ^f + " ^g.<'(^), (4-34) 

Q± = 0. (4.35) 

These values for f2 CT are the only extrema in F a , and a simple check confirms that the upward fluxes are 
maximized at these values. Furthermore, the flux through the upper surface is identical to the downward 
flux for this boundary condition. 

One might assume that the maximal-flux boundary condition is equivalent to applying a radiation 
condition at s = 1. However, this is only true for the kink wave. A radiating upper boundary requires 
that Q a = 0. If the p-mode eigenfunction vanished at the upper surface, the maximal-flux condition would 
indeed be fiy = 0. The fact that the driver is nonzero at the surface ensures that the radiation condition 
and maximal-flux condition are different requirements. 



5. Results 

In the following two subsections we examine several properties obtained from the energy fluxes derived 
in £14.21 In particular, in £15.11 we present the energy flux of tube waves driven into the upper atmosphere by 
the p-mode driving. In £15.21 we examine the damping imposed on the p modes themselves by the excitation 
of tube waves. 
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Fig. 3. — Energy flux of tube waves driven up through the model photosphere summed over all mode orders n < 8. 
The different colors correspond to models with flux tubes with different values of (3: (5 — 0.1 (black), j3 — 1.0 (red), 
and P = 10.0 (blue). 

An examination of the energy- loss rates, equations (|4.26p ~ (14.29j) . reveal that all can be computed di- 
rectly from the p-mode eigenfunctions Q n (w), the eigenvalues n n (v), the mode-energy integral Af n , and the 
interaction integrals In and X±. We compute all of these quantities using a shooting technique with adaptive- 
step-size fifth-order Runge-Kutta integration. The cigcnproblem is solved iteratively at each dimensionlcss 
frequency v and for all mode orders n < 8. Subsequently, once the cigenfunction and eigenvalue have been 
established, all integrals (Af n , I|| , and X± ) are computed during the final iteration of the numerical integrator. 



5.1. Atmospheric Fluxes of Tube Waves 

In ^4. 21 we derived, for a single flux tube, the energy flux of tube waves that passes through the upper 
surface into the upper atmosphere. Since we are interested in comparing the derived fluxes with energy 
fluxes observed within the chromosphere or corona, the energy flux produced by a single flux tube is not 
of particular relevance. Instead, we need the energy flux produced collectively by all of the flux tubes that 
might contribute to the observations. It is easiest and most direct if we examine the aggregate behavior 
of the entire Sun. The total rate of energy deposition £^. ot is simply the deposition rate for a single tube 
E a multiplied by the number of tubes on the solar surface N. If we assume that all of the flux tubes that 
pierce the photosphere arc identical (with the same value of 0, magnetic flux 0, etc.) then the number of 
flux tubes depends only on the surface cross-sectional area of a single tube A s and the filling factor /. Using 
N = ATrR%f/A s and equations (OSjH|OSJ|, one finds 



E\ 



1 (3 Tv 2 gp ujfRl 



\A n \ 



2 + 7 /3(m + l)(,5+l) zl 



K 2 gpoufR 2 Q \A n \ 
(m + l)(/3 + l) zl 



(5.1) 
(5.2) 
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Fig. 4. — Damping rate of / and p modes caused by the excitation of both sausage and kink waves on thin magnetic 
flux tubes. The tubes have f3 = O.f. Two different photospheric boundary conditions are applied: stress free and 
maximal flux. The nulls in the damping rate seen between 5 and 6 mHz for the kink wave under the stress-free 
boundary condition arise because the downward-propagating wave and the wave that reflects off the upper surface 
destructively interfere. 

The final step required to convert these equations into a physically meaningful graph is specification of 
the amplitudes A n of the incident p modes. We use measured p-mode energies to accomplish this. Figure 1] 
shows p-mode energies for harmonic degrees between 4 and 150 as measured bv lKomm. Howe fc Hilll (|2003r ). 



The energies are largely a function of frequency alone and reach a peak value of roughly 2.5 x 10 ergs at a 
frequency of 3.2 mHz. To either side of this peak the energy falls off exponentially with frequency. In order 
to extrapolate to frequencies outside the observed range, we fit these energies with a curve of the form 



Entiv) = aexp - 



(5.3) 
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In this equation, rj is the cyclic p-mode frequency, a is an amplitude, 770 is the frequency of peak energy, d 
is a wing decay rate, and S is a core width. Using a maximum-likelihood procedure we find a = 3.71 x 10 28 
ergs, 770 = 3.16 mHz, S = 0.152 mHz, and d — 0.372 mHz. This fit is shown in Figure [1] with the red curve. 
With this fit we can estimate the p-mode amplitudes A n by equating the mode energies in equation (|3.13|) 
with the fit, 

\A I 2 - H^_fl ^Bt_ (54) 

With this final ingredient we can evaluate the energy-deposition rates appearing in equations (|5.ip and 
(|5.2[) . These rates represent the amount of energy that is driven into the upper atmosphere in the form of 
sausage and kink waves. Figure [2] presents the resulting energy-deposition rates for (3 = 1.0. Figure [3] shows 
the rate summed over all p modes for three different values of (3, representing weak {(3 = 10), intermediate 
(f3 = 1.0), and strong magnetic field ((3 — 0.1). A value of /3 = 1.0 is appropriate for small magnetic elements 
in the photosphere, corresponding to a surface magnetic field strength of 1.2 kG. The other values of (3 
correspond to surface magnetic field strengths of 0.53 kG for /3 = 10 and 1.6 kG for (3 = 0.1. We have 
included the results for all three values of (3 to demonstrate the behavior as (3 changes. In particular, note 
that as the surface field strength decreases (/3 increases) the energy-transfer rate for kink waves decreases, 
whereas the rate for sausage waves increases. For all values of (3 the falloff at both high and low frequency 
is the result of the decrease in the measured p-mode energy spectrum (and hence amplitude A n ) away from 
the peak at 3.2 mHz. 



5.2. p-Mode Damping Rates 

Since the tube waves carry away energy from the acoustic cavity, this energy must come from the p 
modes themselves, thereby damping the solar acoustic modes. We define a damping rate separately for the 
excitation of sausage and kink waves, 

r --s|- (5 - 5) 

Using the energy-deposition rates provided in equations J4.26[) - (|4.29p . we compute the damping rates, 



r || 

LO 



J__L 

10 



A s A" 



fin 



In 



Oi 



S 



A{(3+ l)e|| 4ttR% u 4 



A A? 



fix — 1* 



2 7 (/3+ 1) 4ttR 2 q v 4 



(5.6) 



(5.7) 



In each of the expressions above, the first term in the numerator of the final fraction corresponds to energy 
lost down the tube, while the remaining terms in the numerator arise from energy lost through the upper 
surface. As stated previously, for the stress-free boundary condition the downward term is identically zero, 
and for the maximal-flux boundary condition the upward and downward terms are equal. 
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Figures [4H6] present these damping rates for three values of (3. The most unusual feature within these 
graphs occurs for = 0.1 in Figured] At a frequency around 5.5 mHz, the damping rate arising from the 
excitation of kink waves with the stress-free boundary has nulls where the damping drops to zero. Such 
nulls are to be expected. First, as previously noted, the upward direct wave and the downward direct wave 
generated by the driver possess the same energy flux. Second, the stress-free boundary reflects all upward- 
propagating wave energy back into the convection-zone model. These two facts lead one to the conclusion 
that the reflected wave and the downward direct wave have the same amplitude. There will be frequencies 
at which the phases of these two components are such that complete destructive interference occurs. At such 
a frequency, the da mping rate will vanish b ecause energy is driven neither up nor down the tube. A similar 
result was found bv lCrouch fc Callvi (|l999h . 



6. Discussion 



We have computed the energy flux of sausage and kink waves that are generated in the solar convection 
zone by the buffeting of magnetic fibrils by ambient / and p modes. We treat the fibrils as vertically-aligned 
thin flux tubes. The tubes act as waveguides, much like a fiber-optic cable that ducts the propagation of 
light. Waves are generated locally within the p-mode cavity and freely propagate both up and down the 
tubes — one can easily verify that there are no turning points for our model's stratification. Those waves 
that propagate downward extract energy from the incident p mode and disappear into the convection zone. 
Those waves that propagate upward also damp the p mode and travel through the photosphere along the 
flux tube into the upper atmosphere. Such waves may well be the source of both coronal-loop oscillations 
and the ubiquitous waves observed throughout the corona. 



6.1. Energy Flux of Tube Waves in the Upper Atmosphere 

To determine if p-mode buffeting of small magnetic elements is a viable mechanism to explain the 
waves observed in the corona, we directly compare our derived wave fluxes with the observed fluxes. First, 
we comment that the spectral behavior of our derived energy fluxes is correct by construction. As stated 
previously, the shape of observed coronal wave spectra mimics spectra of p-mode energies — a broad peak at 
3.5 mHz without resonances or nulls. Since our energy fluxes are directly proportional to the p-mode energy 
spectra, we clearly match this observed property well. Now, to determine whether the predicted amount of 
wave energy is reasonable. Table 1 provides the results of summing the tube-wave fluxes generated by all 
p modes with a frequency below 5 mHz. In different columns, we include both the energy flux (ergs cm~ 2 
s _1 ) and energy-transfer rate (ergs s _1 ) for the entire solar surface. 

Table 1: INTEGRATED ATMOSPHERIC ENERGY FLUX 



P E \\/f 

{ergs s -1 } 



0.10 3.81 x lO^ 7 " 
1.00 3.54 x 10 28 
10.0 1.84 x 10 29 



E l} /(4nR%f) 
{ergs cm -2 s -1 } 



E±/f 
{ergs s — 1 } 



6.25 X 10* 
5.82 x 10 5 
3.02 x 10 6 



1.01 x 10^~ 
5.41 x 10 28 
1.35 x 10 28 



E ± /(A7TR%f) 
{ergs cm -2 s — 1 } 



1.65 X 10 6 
8.88 x 10 5 
2.22 X 10 5 



The energy fluxes presented in Table 1 are quite significant. |Pe Moortel &: Rosnerl (|2007j ) estimate that, 
depending on the mass density of coronal loops, the energy carried by loop oscillations is E « f c x 3 x (10 25 — 
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Fig. 5. — Same as Fig. |4j except that the thin flux tubes have (3 = 1.0. 



10 26 ) ergs s 1 , where f c is the coronal filling factor. iTomczvk et alj (|2007f ) estimate that the upward energy- 



flux of waves in the high corona is 10 ergs cm 
area. 



or E « 10 ergs s if multiplied by the solar surface 



If we adopt a typical value of f3 = 1.0 for magnetic elements in the photosphere and a value of / = 0.01 
for the photospheric filling factor, we find that our derived energy-transfer rates are 3.5 x 10 26 ergs s _1 and 
5.4 x 10 26 e rgs s -1 for the sausag e and kink waves, respectively. While these numbers s lightly exceed the 
estimates of lDe Moortel fc Rosnerl ([20071 ) and exceed the estimate of lTomczvk et al.1 (|2007l ) by several orders 
of magnitude, we must keep in mind that our theoretically derived fluxes are the maximum possible flux 
that can be driven upward at the photospheric level. Reflection and absorption within the chromosphere 
and transition region are ignored. Clearly sufficient energy can be converted from p modes to tube waves 
to explain the observations; however, whether the attenuation of such waves allows them to reach high into 
the corona remains to be seen. We remind the reader that flux tubes can act as waveguides. The acoustic 
cutoff frequency of the nonmagnetized atmosphere is not the relevant quantity to determine whether the 
tube waves are propagating. Instead, the sausage and the kink waves have separate cutoff frequencies which 
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Fig. 6. — Same as Fig. |4j except that the thin flux tubes have (3 = 10.0. 



can be lower than the acoustic cutoff of the surrounding m edia, depending on the properties of the tube 



(sec 



Roberts fc Webblll978t iMusielak fe Ulmschneiderll200lh . Therefore, inclination of the field may not be 



required for magnetosonic waves to "tunnel" through the photosphere and temperature-minimum region. 



6.2. Damping and p-Mode Line Widths 

The excitation of tube waves extracts energy from the incident p mode that drives the oscillations. 
Figures HHB] show the resulting damping rate. In order to determine whether these damping rates are 
significant, in Figure[7|we plot t he damping rates for = 1 .0 and a filling factor of 0.01. We then overlay the 
p-mode full widths measured by Komm. Howe fe Hilll ( 2003h . From this figure, one can see that the excitation 
of tube waves is likely to be a significant damping source only for the lowest order p modes. The damping 
is also a strong function of the boundary condition. The stress-free boundary produces significant damping 
only for the / mode, and most of the energy is converted into kink waves. For the maximal-flux boundary, 
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Fig. 7. — Same as Fig. [4] except that the value of (3 is 1.0 and measured full widths of the p modes have been 
overplotted with diamonds. Each mode order is plotted with a different color, using the same color s chem e as in the 
theoretical curves: / (black), pi (red), etc. The full widths were measured by Komm. Howe fc H ill (2003) and only 
orders 1-7 are shown. 



markedly more energy is converted into sausage waves than kink waves, and the damping is important for 
n < 5. 



6.3. Conclusions 



We have demonstrated that p-mode buffeting of thin flux tubes can easily generate the flux and spectral 
dependence of waves observed in the corona, in the form of either longitudinal sausage waves or transverse 
kink waves. However, we only show that these large fluxes are possible at the photospheric level. Therefore, 
the remaining problem is one of transmission: Can the tube waves successfully propagate high into the corona 
without being absorbed or reflected at lower levels in the atmosphere? We plan in future efforts to include 
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a realistic model of the temperature minimum and low chromosphere, thereby accounting for the flaring of 
flux tubes with height and their eventual merger into a magnetic canopy. However, the semianalytic method 
used here is likely to prove of limited usefulness once the model atmosphere and model flux tubes become 
sufficiently complicated. Numerical simulations of the interaction between p modes and flux tubes will 
probably provide the ultimate answer. Such models have been previously built to investigate the absorption 
of acoustic waves by sunspots; however, these models have by explicit construction eliminated upward 
propagation of the resulting tube waves. The inclusion of a chromosphere and realistic upper boundary 
conditions is needed. 
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